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MACROSCOPIC NON-UNIQUENESS AND LIMITS OF 
HAMILTONIAN DYNAMICS 


S. DOSTOGLOU AND JIANFEI XUE 

Abstract. We construct explicit examples of spontaneous energy generation 
and non-uniqueness for the compressible Euler system, with and without pres¬ 
sure, by taking limits of Hamiltonian dynamics as the number of molecules 
increases to infinity. The examples come from rescalings of well-posed, de¬ 
terministic systems of molecules that either collide elastically or interact via 
singular pair potentials. 
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1. Introduction 

Non-uniqueness for weak solutions of hydrodynamic equations is well known. 
Examples include the construction of V. Scheffer [Sch] and A. Shnirelman [Sh] 
of non-trivial weak solutions of (incompressible, two-dimensional) Euler equa¬ 
tions with compact time and space support, and the work by C.DeLellis and 
L. Szekelyhidi [dLS] showing that non-uniqueness (of the incompressible and com¬ 
pressible Euler equations in dimension greater or equal to two) persists even under 
“admissibility” conditions. [D] is a standard reference on the non-uniqueness of 
weak solutions of hyperbolic conservation laws in general. 

In an attempt to investigate the origin of this behavior, we adopt here the 
point of view that hydrodynamic equations are the result of averaging microscopic 
evolution equations (cf. [M], p.81, and [B], Part I, §20) to construct explicit ex¬ 
amples of spontaneous macroscopic energy generation and non-uniqueness for the 
compressible Euler system, with and without pressure, as limits of Hamiltonian 
dynamics. Our examples are rescaled limits of well-posed, deterministic systems 
of molecules that either collide elastically or interact via rescaled, singular pair 
interaction potentials, at the limit of infinitely many molecules, cf. C.B. Morrey’s 
work [Mor]. For each moment t and finite N, the positions and velocities of the 
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1 ^ 

molecules define the probability measure dv) := — E <^(xfc,Vfc)(dx,dv). 


k=l 


In all examples here, converges weakly to Mt as — )■ oo and for each {t, x) 

the macroscopic density is given by the first marginal of Mt and the macroscopic 
velocity by the barycentric projection of Mt at x with respect to this marginal. 

The first part of this article, consisting of Sections 3 and 4, is centered on an ex¬ 
ample showing spontaneous generation of macroscopic velocitiy. The microscopic 
systems start with groups of motionless molecules and a single molecule, macro- 
scopically undetectable, initially at a sufficiently large distance from the group, 
moving towards the group. Macroscopically, the limit of these flows describes 
a line segment in at rest for t G (—oo,0], which splits into two equal parts 
moving away from each other with velocities ±1 as soon as t becomes positive. 
The macroscopic velocity and the macroscopic density from Mt turn out to be a 
weak solution of the 2-dimensional presureless Euler for all t in R. This solution 
is macroscopically as “inadmissible” as those of Scheffer and Shnirelman in that 
kinetic energy is spontaneously created at to- (Microscopically, total energy is, 
of course, conserved.) As Hamiltonian flows are time reversible, in Section 4 the 
flows M^^'^ are reversed to produce a solution to the 2-dimensional presureless 
Euler that does decrease energy. When compared to an elementary transverse 
flow, this provides an example of non-uniqueness for the presureless Euler under 
the admissibility condition of non-increasing energy. 

In the second part of this article. Section 5 provides an interpretation, via a 
microscopic derivation, of the well known non-uniqueness of the Cauchy problem 
for the 1-dimensional Euler system. We show how three moment equations derived 
from the transport equation 


( 1 . 1 ) 


dtMt + vdxMt = 0 


can result in the 1-dimensional Euler system. The main point here is that two 
flows of probability measures solving the same transport equation, even if their 
moments coincide at t = 0, in general will not have identical moments for all 



4 


S. DOSTOGLOU AND JIANFEI XUE 


later times. Indeed, we construct two limit measures Mt and Mt both solving 
(1.1) and resulting in the 1-dimensional Euler system. At t = 0, both Mt and Mt 
give the same macroscopic density, velocity, and pressure. Macroscopically, the 
solutions produced by Mt and Mt can be pictured as a segment of two and three 
layers, respectively, on top of each other moving freely, see Figures 11 and 13. The 
solutions in Section 5 are surrounded by vacuum (zero density). 


2. Preliminaries and Notation 


2.1. Measure theory. Recall that a sequence of finite measures Mn{dx) con¬ 
verges weakly to a finite measure M{dx) if for any f{x) continuous and bounded 
J f{x)Mn{dx) J f{x)M{dx), n ^ oo. We then write Mn M. 


For f : X ^ Y measurable and M a probability measure on X the push-forward 
measure f#M of / on F (the distribution measure of the random variable /) is 
{f^M){B) = M{f~^{B)). We often write fM for this push-forward. 


If M is on its first marginal will be (7ri)^M, for tti : —)■ M“,7ri(x, v) = x 


Tfx('^v)/r((ix) is a shorthand for the measure / I—)■ / ( / /(x,v)Mx((iv) ) ^(dx) 


The disintegration of M{dx,dv) with respect to its first marginal ;u(dx) is the 
unique, up to a /U-measure 0, family Mx(dv) such that M (dx, dv) = f M^{dv)fi(dx.). 


The barycentric projection of this disintegration is v(x) = J vMx(dv) for x in 
the support of /i, v = 0 otherwise. For details see [AGS], Section 5.3, or [DJX], 
Section 3.1. 


2.2. Finite systems. A system of N molecules in will be described by the 
positions and velocities of the molecules, (xfc(t), Ufc(t)), 1 < A: < A^, evolving via 
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Hamiltonian dynamics with pairwise interaction <ho-(p) of finite range a: 


( 2 . 1 ) 


d , , 


dt 


^k{t) 


Ufe(t), 


1 

N 


N 


i=i 


Xj(t)l) 


Xfc(t) 

Xfc(t) 


Xj(t) 

xi(t)r 


Following Morrey [Mor], we shall take <ho-(p) = ‘h ^ for some <h ; (0, oo) —)• [0, oo) 
satisfying: 


(2.2) lim ^>(r) =+ 00 , < 0, $" > 0, $(r) / 00 < r < 1. 

r—>-0 

For each N, suppose that a system ^x^^^(t), u^^^(t)^, k = l,...,iV evolves 
according to (2.1). Of central importance will be the corresponding t-family of 
probability measures on 

1 ^ 

(2.3) Mi’^\dy.,dw) := ^^^Md:i^,dyr), f > 0, or f G R. 

fc=i ' ' 

When converges weakly to some Mt, it is crucial to note that the empirical 

measure formed by neglecting a single molecule converges weakly to the same Mf. 
(In fact, neglecting o{N) number of molecules has the same effect.) In this sense, 
any single molecule is macroscopically invisible. The construction in Section 3 
relies heavily on this observation. 


3. Spontaneous Macroscopic Velocity Generation from Hamiltonian 

Dynamics 

This section presents an example of a microscopic Hamiltonian flow with macro¬ 
scopic limit, as N ^ oo, that shows spontaneous velocity generation. The micro¬ 
scopic systems start with groups of motionless molecules and a single molecule, 
initially at a sufficiently large distance from the group, moving towards the group 
with large velocity. For t < 0, as V —)■ oo, the moving molecule is invisible and the 
macroscopic system is motionless. However, as the moving molecule starts inter¬ 
acting with the group at t = 0, its energy is transferred to the rest of the system in 
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such a way that all other molecules acquire speed 1 to create macroscopic velocity 
for t > 0. 

There are similarities here with Lanford [L], pp. 50-53, although Lanford works 
with an infinite system of hard balls that always remains discrete, rather than 
the limit of finite Hamiltonian systems with interaction, and he does not obtain 
hydrodynamic equations. 

Throughout this section we use Qt for the segment 

(3.1) {{x, y) : 0 < X < 1, y = t} C 


and At{dx) for the normalized 1-dimensional Lebesgue measure on Qt. 


Theorem 3.1. For each N £ N, there exists > 0 and (t), 
k = 1,..., N, solution of the Hamiltonian system (2.1) with interaetion for 
all t G M, sueh that for all t G M, the sequence of empirical measures 

(3.2) Mf)(dx,dv) := 

k=l ' ' 

eonverges weakly, as N ^ oo, to 


(3.3) Mt{dx,dv) = 


Ao(dx) (g) (5(0,0) (d'v) t <0 

^-At((ix) (g) (5(0,1) (dv) l-A_t{dx) (g) (5(o,-i) (dv) t > 0. 


The proof of this theorem occupies the rest of this section. For the moment, 
note that the first marginal (macroscopic density) of Mt(dx, dv) in (3.3) is 

t < 0 


(3.4) 

If we disintegrate 

(3.5) 


Ao(dx 

Mt(dx) = { I 


1 


-At(dx) -A_t(dx) t > 0, 


Mt{dy.,dv)= / Mt^^{dw) yt{d: 


X 


then 

(3.6) Mt^^{dv) 


TQo(x)d(o,o)(dv) t<0 

XQt (x)d(o,i) (dv) XQ-t (x)d(o,-i) (dv) t > 0. 
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t < 0 


t > 0 


Figure 1. Macroscopic flow of Mt in Theorem 3.1. 


Notice that in (3.6) we only needed to specify for x in the support of 

Ht^dx). The macroscopic velocity is the barycentric projection of this disintegra¬ 
tion: 


(3.7) u(t,x) 



(0,0) t<0 

XQt(x) • (0,1) +yQ_,(x) • (0,-1) t>0. 


The macroscopic density (3.4) and velocity (3.7) show clearly a macroscopic ve¬ 
locity generation (see Figure 1): before t = 0, the macroscopic system stays at 
rest, while, starting at t = 0, two equal mass fronts split and move away from 
each other with velocity ±1. The sudden increase of macroscopic kinetic energy, 
of course, comes from interaction with an invisible molecule as we will see in the 
proof of Theorem 3.1 (subsections 3.1, 3.2, and 3.3). In subsection 3.4 we examine 
the macroscopic hydrodynamic equation solved by the density (3.4) and velocity 
(3.7). 


3.1. Interaction with one particle at rest. Start with two identical molecules 
P, Q interacting with potential as in (2.1). Denote the positions and velocities 
of P, Q as xp = {xp,yp), xq = {xQ,yQ), vp, and vq. Consulting Figure 2, let D 
be the disc with center (xq, yo) and radius r > 0 and assume that at t = 0 

(1) (xp, yp) G D and xq = xq + d with d > r + a, i.e. P is inside D and Q is 
on the vertical line x = xq + d. 

TT TT 

(2) Vp = u(cos(/>, sin(/)) with —— < cj) < —, v > 0 and vq = (0,0), i.e. P 
moves with speed v and Q is at rest. 
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Figure 2. The initial disc D and the segment I for (j) < 0. 


We say that there is interaction between P and Q whenever their distance is 
smaller than a. Since Q is at rest at t = 0, there will be no interaction between P 
and Q as long as P is inside D. The following lemma on the interaction between 
P and Q is the building block of the rest of this section. 


Lemma 3.2. Let P, Q he as above: 

(1) For any 9 in ^ j there exists yq such that P and Q will eventually in¬ 
teract (i.e. P and Q will interact at some time t>t)), and after interaction 
P and Q will move in directions perpendicular to each other with constant 
velocities Vp = u cos 0 (cos 0', sin(^') and Vq = u sin0 (sin(^', — cos 0'), re¬ 
spectively, where (f)' = 4> + 9. 

(2) If interaction takes place then yq satisfies 


VQ 


(yo + d tan 0) | < 


r a 
cos (p 


( 3 ) 


Whenever P and Q interact, they are both inside the disc with center (xq + 

r + fj 

a, yo + d tan (p) and radius -h 5a. 

cos (p 






MACROSCOPIC NON-UNIQUENESS 


9 



Proof. (1) Consulting Figure 3 (which is [LL]’s Figure 17, p.47, in our notation), 

for 6 the deflection angle from vp to Vp, conservation of momentum and energy 

/ TT 7r\ 

gives the formulas of Vp and Vg. That any 0 in is attained by some 

Uq follows from Corollary A.3 in the Appendix and the formulas in [LL], §13 that 
show how to transform from motion in a central field to a system of two molecules. 

(2) Let S be the strip between the two lines tangent to D and parallel to vp, 

Sa the set of all points with distance smaller than a from S', and I be the interval 

of intersection of with the line x = xq + d, see Figure 2. Then if Q has second 

coordinate anywhere out of /, P ignores it and continues with unaltered velocity 

Vp. Elementary geometry shows that I has midpoint yo + dtancf and half-length 
r -|- cj 

--, consult Figure 2. 

cos (j) 

(3) By Lemma A.l in the Appendix, when P and Q interact, their interaction 

4cj 

time is less than — and by conservation of energy is positive) the speed of 

V 

Q will never be more than v during interaction. Therefore, during interaction 


Q travels less than Aa, i.e. it stays in the disc centered at (xq + d,yQ + dtancp) 
r a 

with radius - - + Aa. As the distance between P and Q is always less than cj 

cos 4> 


during interaction, P is always inside the circle centered at (xq + d,yo + dtancj)) 


with radius 


r + a 

-7 + 5cJ. 

cos (j) 


□ 


3.2. A system of molecules on the plane. We describe now a system consisting 
of A-l-1 molecules P, Qfc,/? = !,...,A where P interacts (only once) with each Qk 
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(in the order of increasing k) and interactions are independent (P does not interact 

with Pj, j ^ k, when interacting with P^, and there is no interaction between the 

Qfe’s). In addition, the moment before interacting with Qk the speed of P will be 

greater than 1 and the speed of Qk after interaction will be 1. 

We use dk for the deflection angle of P due to the interaction with Qk- Assume 

k 

that before interacting with Qi, P moves along the x-axis. Then (j)k = '^^Oj will 

i=i 

be the angle from the x-axis to the direction of the velocity of P right after its 
interaction with Qk- The angle from the x-axis to the direction of the velocity 
of Qk after its interaction with P will be denoted by 4>k- By Figure 3, (pk = 

Lemma 3.3. For N gN fixed and k = 1,2,... ,N, let 

k 

(3.8) Ok = (-1)^ arcsin . </>fc = ^ Oj- 

^ i=i 

Then 

(1) (pk < 0, when k is odd and (pk > 0, when k is even, 

(2) \(pk\ < \(pk+2\, 

(3) \(pi\ = |6»i| < ^ and \(pk\ < \0k\ < for k>l. 

Proof of Lemma 3.3. For (1), observe that the 9k s start negative, increase strictly 
in absolute value and alternate sign. Therefore for k odd and k > 1 

(3.9) (pk = Gi + (02 + Os) + ... + {9k-i + Gk) < 01 < 0, 
whereas for k even 

(3.10) (pk = (01 + 02) + ... + {0k-i + ^fc) > 0. 

For (2), notice that 0k+i + 0k+2 always has the same sign as (pk, hence 


(3.11) 


I0fc+2| — + l^fc+1 + (^k+2\ > \(pk\ • 
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Figure 4. The angles (pk and 4>k- Observe how the even/odd (pk’s 
and the even/odd ^^’s fall into four non-overlapping sectors. 



(3.14) \(j)k\ = (t)k = (j)k-i + 6'fc < 4 < ^- □ 

Lemma 3.3 shows that the even c/i^’s are positive, increasing, and never more 

than 7r/4 (and therefore the even (p^s are negative, increasing, and never more 
than —7r/4), whereas the odd (f>k’s are negative, decreasing, and never less than 

—7r/4 (and therefore the odd (p^s are positive, decreasing, and never less than 

7r/4). Figure 4 summarizes the behavior of (pk and (pk- 
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In the description of the interaction of P and Q^, for (jij as in (3.8) (assuming 
00 = 0 ), the point 


(3.15) 


{xk^Vk) = 


k 1 v-^ 

—, — > tan( 

N' N ^ 


3=0 


will play the same role as (xq + d, yo + dtan0) in Lemma 3.2. The segments and 
half-lines 


(3.16) 

Vk := {{x,y) ■■ Xk<x < Xk+i, y = {x - Xk)tan(pk + Vk} , k = 1, ... ,N - 1, 

Vn ■■= {{x,y) -.XN <X, y = {x- XAr)tan0Ar + yN} , 

Qk ■= \^ix,y) : X > Xk, 2 / = (x - Xfc) tan yfcj , k = l,..., N, 

will be useful in describing the trajectories of P and of each Q^, respectively, see 
Figure 5. Define the distance between any two of these sets as 


(3.17) 


d{A,B) :=inf{||a-6|| : a £ A, b £ B} . 
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Lemma 3.4. Let m,n = 1,2,..., N. For Qm, Vn as above, 

d {Qmi Qn) ^ ■) rn ^ n, 

(3.18) ^ 

d{Qm,Vn) > ^, m<n. 

Proof of Lemma 3.4- Recalling (3.16), we use here “right half plane of Qn^ to 
mean the half-plane to the right of the y-axis defined by: is the positive y-axis 

when (f)n is positive; is the negative y-axis when <fn is negative. 

Observe first that for any fixed n, the point {xn,yn) is always in the right half 
plane of Qm for all m < n : this holds by the relation of the angles cfi to the angles 
4>j, see Figure 4 and Figure 5. 

To get the first estimate in (3.18), it suffices to consider n > m. If n — m is 
even, then the angle of Qn (i.e. (/>„) is of smaller absolute value than the angle of 
Qm (i.e. (f>m)- If n — m is odd, then the angles of Qm and differ by more than 
7r/2. In either case the point on Q„ closest to Qm is {xn,yn)- 

Similarly, the angle of Vn (i-e. (pn), is always of absolute value smaller than the 
angle of any Qm (i.e. (pm)- Therefore the point on Vn closest to Qm, for m < n, is 
{Xm yn)- 

Now it suffices to notice that the distance from (x„, yn) to each Qm is greater or 
equal to \Vm\ which is clearly bigger than — (consult Figure 4 and Figure 5). □ 

We are now ready to establish the evolution ol P,Qi,..., Qn- 

Proposition 3.5. For each iV G N and gn < —:=—^- zntk, consider the 

2\/2 A^(A^ + 3)3/2 ’ 

system P,Qi,... ,Qn with interaction <l>o-jv. Then there exist yq^, ’s, k = 1,..., N, 
such that the system evolves as follows: for t < 0, 

P{t) = tvp, vp{t) = (^N + 1,0^ , 

(3.19) X ^ X 

Qk{t) = y—,yQ^] , VQjt) = (0,0), k = l,...,N, 


and for t > 0 
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(1) There exist times 0 < < t'{ < t '2 < 1 2 < ... < t'j^ < t'^ such that for 

any 1 < A: < P starts to interact with Qk at t = t'j. and completes this 
interaction at t = t'l. 

(2) For any 1 < k < N, the molecule Qk does not interact with any other 
molecule for t <t'^ or t > t'l and its velocity is given by 


(3.20) 


^Qkit) 


/ 

(0,0) t < 4 

(sm|(^fc|,(-l)*'+^cos(/>fc) t>t'l, 


for fk as in (3.8). 

(3) The velocity of P satisfies 


(3.21) vp(t) 


iVWTT, 0) t < 4 

< VN + 1- k (cos (pk, sin (fk) t'l < t < t'i.^.^^,k = 1,2,..., N - I 
(cos (()7v, sin (/)jv) t>t'ff. 


for (fk as in (3.8). 

(4) During the time interval [4)41 far 1 <k < N the molecules P and Qk are 
in the disc of center {xk,yk) as in (3.15) and radius given recursively by 


(3.22) Tfc =---h5(TAr, ro = 0. 

coschk-i 


In partieular, 


(3.23) 


Vk < 2V2{N + 3)^/^crAr. 


Proof of Proposition 3.5. For all t < 0 and any choice of i/Qj,, k = 1,..., take 
^^p{t) = (VA^ + 1,0), P{t) = Up, Qk{t) = (^■^,yQ^'^, and ^rQ^{f) = (0,0). For all 

aN < —, it is clear that P,Qi,..., Qn solve the Hamiltonian system for t < 0 (as 
there is no interaction). We now specify 2 /Qj,’s for the evolution when t > 0. 

Applying Lemma 3.2 for xq = yo = 0, r = 0, <() = 0, u = ^/N + 1, d = 1/N and 
6 = 61 = 4>i = — arcsin(l/\/Al + 1) there is yg^ such that P will interact with Qi 
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Figure 6. Schematics of the system of Proposition 3.5 after three 
interactions. The radii discs are not up to scale. 

and after interaction 

(3.24) vp = \/iV (cos sin ^i), = (—sin cos . 

In this way, the position of Qi, depending on ajsf, is now determined. The whole 
interaction, according to Lemma 3.2, takes place in the disc of radius ri = Gittv 
and center (1/iV, 0). Let [tijt"] be the time interval of this interaction. Preparing 
for the next interaction, make a new choice of ajsf so that ri = GuAr < 1/A^, and 
note that everything in this first step still holds for the new choice of (Tat. 

For induction, fix /c G N and assume that ri, ■ ■ ■ ,rk satisfy (3.23), and therefore 
rj < 1/N, j = 1,... , k, for all aN small enough. Further assume that yq^,. ■ ■ ,yQ^, 
t[,... ,4, t'(,... ,tl, VQ^(t),..., vp{t), for t < t'l, have all been determined 

and satisfy (3.21) and (3.20). 

Apply Lemma 3.2 for (xo,yo) = ixk,yk), for {xk,yk) as in (3.15), r = Vk, 
4> = (l>ky V = \/N + 1 — k, d = 1/N and 6 = 6 k+i as in (3.8), to find that Vk+i is 
determined by formula (3.23), to determine the times and the 

velocities vp{t), VQj,_|_^(t) for t G that will satisfy (3.21) and (3.20). 

Therefore Qk+i is always in the r^+i-neighborhood of Qk+i-, as defined in (3.16). 
Choose cjat so that is smaller than 1/N. Using Lemma 3.4, Qk+i does not 
interact with Qi,... ,Qk during the interval (—oo, 
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For (3.23), rewrite first (3.22) as 


Ffc — sec (f>k—i'>’k—i + sec (pk—i^N + Scttv 
(3.25) k—l k—l k—1 /k—1 k—1 

= sec (pjro + ^ n sec (pmC^N + ^ sec + 1 j Sctat 

j=0 j=0m=j yj=lm=j 

and, using ro = 0 and \(j)j\ < \6j\ (Lemma 3.3), estimate this by 


k—1 k—1 


k—1 k—1 


(3.26) sEn sec + En sec 0m + 1 j 5iTAr, 

j=0m=j \j=lm=j 


and then, increasing A: to and using (3.8), estimate the same by 


N-l 

sE 

j=0 


VN+2-j 

V2 


' N-l 


-aN 


+ E 


VN + 2-j 


,i=i 


^/2 


+ 1 5cJAr 


(3.27) 


N-l 


N+2 


< 


< 


j=0 V j=3 

/*7V+3 

3\/2cJAr J \/xdx < 2V2{N + 3)^^‘^aN. 


In particular, cjat < — j =—implies < 1/N for all k. 
2v2 N[N + 3)'*/^ 


□ 


Remark 3.6. Notice that, for each N, Proposition 3.5 provides examples of the 
general theory of Gal’perin and Vaserstein, [G] and [V], according to which, for fi¬ 
nite range interactions, molecules evolve by eventually separating into independent 
clusters. Each cluster here consists of a single molecule. 


3.3. The limit system as iV — )■ oo. In the notation of Proposition 3.5, let T)(r := 
N 

{t" — t'j), the time during which P interacts with some Qk- Then Tfj = t'fi — 
is the time during [0, t'fi] when P is not interacting at all. 


Proposition 3.7. t'fi ^ 0, as N ^ oo. 
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Proof. According to (3.21), the speed of P at is VN + 2 — k. Then, by Lemma 
A.l, 


(3.28) 


N 

U/<E 


4(TAr 


^ VN+ 2-k 


N+l 
4<TAr ^ 


k=2 


1 

Vk 


After the interaction of P with Qfc is complete, P moves with speed VN + 1 — k, 
forming angle (fk with the x-axis. The distance that P will travel until its 
interaction with Qk+i begins, satisfies 


(3.29) 


dk < 


1 1 

< 


N cos ifk N cos 9 k 


TT 


cf. Figure 2. Recalling that \9k\ < — from Lemma 3.3 gives 


7V-1 


T’;f< 


T —^ 

.4—/ A7 r.n 


(3.30) 


k=0 


N cos Ok VN + l-k 


■ N-l 


<^y 

- N ^ 


y/2 1 

~N 2^ TTf.' 


7V-I-1 


^ ^ VN + l-k A' Vfe' 


This and (3.28) imply 
(3.31) 


tjv < 


Af-l-l 




Af-l-1 


1 

k=2 

as N ^ oo. 


As -y= < 2V A" + 1, and for as in Proposition 3.5, we conclude that t'ff —)■ 0 

□ 


Proposition 3.8. max^yQi^ — )• 0 as iV —)■ oo. 

Proof. Noting that yq^, is the second coordinate of Qk before t = t'^, whereas yk 
is the second coordinate of the center of the fc-interaction disc, it follows from the 
definition of and (3.23) that 


(3.32) 


IVQkl < \yk\ + ‘2V2{N + 3)^^^CJAr. 
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For the second term on the right use ajsf as in Proposition 3.5 and estimate the 
first term as 


(3.33) 


as N ^ oo. 


\yk\ = ^ 


k-l 


tan( 


m=0 

k-l 


k-l 


- jV ^ tan|(/.^ 


m=0 

k-l 


^ ^ / 


N 

k- 

E 


m=0 
k-l 


m=0 

N-1 


arcsin , 

V VN + 2 - m 


< 


E 


N^N + l-m ^giVViV + l-m 
1 

— > ^ ^ 0, 

N ^ y/m 

m=2 ^ 


□ 


For each fixed N, writing v = {vx,Vy) and following (2.3), set for t G 


(3.34) 

Mj (dx, dy, dvx, dvy) 

1 ( ^ 

= ^ 1 I S{P{t),vp{t)){dx, dy, dvx, dvy) + E '^(Qfc W.VQ;, W) ^^2/) 

^ V k=l 

The crucial observation in the following proposition is that, due to the factor 1/iV, 
no single molecule shows as N ^ oo, but its interaction with many other molecules, 
if their number is of order N, shows macroscopically. 

Proposition 3.9. As N ^ oo, and for as in Proposition 3.5: for t <Q, 

(3.35) M^^^^\dx,dy,dvx,dvy) ^ X[o,i]{x)dx ^ 5o{dy) (g) 5(Qfl){dvx,dvy), 



and for t > 0, 

(3.36) 

Mj {dx, dy, dvx, dvy) 

^St{dy) (g) 5{o^i){dvx,dvy) + ]^5-t{dy) (g) 5{Q_i){dvx,dvy) 


^ Xlo,i]ix)dx 
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Proof. It suffices to check the statement on the integrals of bounded Lipschitz 
functions, see [AGS], page 109. For this, for / : x —)■ M bounded and 

Lipschitz 


(3.37) 


/ {x,y,Vx,Vy) Mj:^^^\dx,dy,dvx,dvy) 

1 1 ^ 

= WTT^ vp(t)) + f VQ,(t)). 


Since / is bounded, the first term vanishes as N ^ oo. The rest of the proof 
examines the convergence of the second term. 

Fix any t < 0. Recalling (3.19), 

1 Af ^ ^ / k \ 

k=l k=l ^ ^ 

For Lf he the Lipschitz constant of /, and using Proposition 3.8, 


(3.39) 


^ ^ / If \ 

FTT^'^ 7 ~ N + i 

k=l ^ ^ k=l ^ 


< - Lf max \voA —^ 0. 


By the definition of the Riemann integral, 


(3.40) 


Therefore 




/ (x, 0,0,0) dx. 


(3.41) f f {x,y,Vx,Vy) Mj;^~^^\dx,dy,dvx,dvy) -f f f {x, 0,0,0) dx. 

4r4 Jo 


This is exactly (3.35). Now fix t > 0. By Proposition 3.7 there exists A^i such that 
for all N > Nf, t'ff < t, i.e. for each time we can choose N large enough so that all 
interactions have already happened and all molecules are moving at time t, and 
are moving with their terminal velocities. We consider such Al’s only. According 
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to Proposition 3.5, and since now t > t'^, 


(3.42) 

For ctat = 
OtN 


XQkii) = XQk{4) + {t- 4) VQ^^a:{4), 

yQk4) = yQki4) + 4 - 4) VQk,y{4)- 

N — '/n] , the integer part of iV — \/]V, and by (3.20), for any 1 < k < 


'>^Qk,x{4)\ =sin|())fc| < sin|6»fc| < 


\/N + 2 — ajsf ’ 


(3.43) VQ^^y{4) - (-1) ^ = |cos(/>fc - 1| < |sin())fc| < 


y/N + 2 — Q!7V ' 


XqM) - 


<rk, \yQki4)\ <\yk\+rk- 


Therefore for 1 < A: < a^, by (3.23), Proposition 3.7, and Proposition 3.8, 

k 


XQkit) - 


N 


< 


XQki4) - 

<rk + 


A 

N 

t 


+ 4~4) 4x,Qk 

0 , 


(3^44) \/iV4-2^-7rjv 

yQk4) - (-l)^+H < \yQ44)\ + t VQ^,y{4) - (-1)^+^ + 4\vQu,y{4)\ 

t 


< \yk\ +rk + 


y/N -\- 2 — aj\j- 

Since / is Lipschitz, (3.43) and (3.44) imply that 
(3.45) 


+ 4 —^ 0 - 


1 


N + . 

k=l 

For Cf = max I/I 
1 


On oiv / T \ 

Y E / «<=(*).VO.W) - E / X' 

k=l k=l ^ ^ 


0 . 


(3.46) 


N 


N + 1 
1 


f {Qk{t),^rQ^{t)) 


fc=ajv+l 
N 


^ ^ N -ON ^ ^ 

^ Ar , 1-^ 0, 


A^ + 1 


N + 1 


^ /(^A,(-i)^+it,o,(-i)^+i) 


k=aN-\-l 


<Cf- 


N — a/v" 
A^ + 1 


0 , 
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therefore, 

(3.47) 


1 


N 


N 


N + 


■Y'^f{Qk{t),VQ^{t)) 


k=l 


fc=i ^ 


0 . 


By the definition of the Riemann integral, 
(3.48) 


N 


] ^ / k \ ^ 


which implies (3.36). 


□ 


With 


(3.49) 


= Qk{t), =\rQ^{t), k = l,...,N 


iN+l)u\- 


^^N+l\t) = P{t), U^N+l\t) = vp(t). 


> —OD • 

/*oo 


Theorem 3.1 follows immediately from Propositions 3.5 and 3.9. 

3.4. Macroscopic equations. We now examine the hydrodynamic equations for 
Mj(dx, dv) as in Theorem 3.1. It is easy to check that for any x) G C'“(M x M^) 

(3.50) 

/ OO p 

/ Vx(/>(t, x) • vMt{dx, dv)dt = 0, 

-OO Jw^ 

/ OO /* 

/ Vx0(t,x) • vvMt(dx,dv)dt = 0. 

-OO 

Using disintegration (3.5), for ^t(dx) and u(t,x) as in (3.4) and (3.7), we rewrite 
(3.50) as 

/ OO p poo p 

/ dt4’{t,x)iJ,t{dx)dt+ / Vx(/>(t,x) • u(t,x)^((dx)dt = 0, 

-OO J — OO ^R^ 

/ OO ^ 

/ dt4’{t,x)\i{t,x)iJ,t{dx)dt 

-OO t/R2 


+ 


Vx(/>(t,x) 


V (g) vMt^:x.{dv) I fj,t{dx)dt = 0. 
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Notice that at each t,x the Mt^y^{dv) is singular, therefore 


(3.52) 

Then (3.51) becomes 

/ OO P 

/ dt(p{t,x)fit{dx)dt + 

■OO 


V (g) vMt^x{dv) = u (g) u. 


Vx</>(t,x) • uij,t{dx)dt = 0, 


(3.53) 


POO P POO P 

/ / dt(p{t,x.)unt{dx)dt + / / V:x_(p{t,x.) ■ uu^tidx)dt = 0. 

J—OO J—ooJm? 


In other words (fitidx), u(t, x)), t G M solves weakly two dimensional Euler system 
without pressure: 

dtUt + div{u fit) = 0, 

dt{ufit) + div(u 0 u fit) = 0. 

For the naturalness of measure solutions in the presureless Euler system, see [ERS], 
p.354. 


(3.54) 


Remark 3.10. The trivial solution Jitidx.) = Ao((ix),u = (0,0) also solves (3.53) 
for all t, and coincides with {fit,u) for t < 0. Note that (^t,u) is not “energy 
admissible” sinee the kinetic energy of {fit, u) inereases in time: 

f O f O fot<0 

(3.55) / \u\^fit{dx) = / |vpMt((ix, dv) = < 

4r2 Jr* I 1 t > 0. 

A solution to (3.53) with deereasing energy ean be obtained by reversing the diree- 
tion of time, as in the next seetion. The value of the eonstruetion in this seetion lies 
in the microscopie, Hamiltonian interpretation of spontaneous veloeity generation 
in weak solutions of hydrodynamic equations as in [Sch], [Sh]. 


4. Time Reversal and Macroscopic Non-Uniqueness 

4.1. Reverse flow with decreasing energy. We now reverse time in the con¬ 
struction of the previous section to establish macroscopic non-uniqueness in the 
class of energy decreasing solutions. It is standard that for ^x^^^(t), u^^^(t)^ a 
Hamiltonian flow as in Theorem 3.1 the reverse flow (x.^j^\—t),—u^f^\—t)'\ also 
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t < 0 t>0 

Figure 7. Macroscopic flow of (4.1). 


solves the Hamiltonian system (2.1). Roughly speaking, for each N the reverse 
system consists of N molecules moving with speed 1 for t < 0. At t = 0, through 
interaction, one of the N molecules gathers all the energy from the rest N — 1 
molecules and leaves the rest of the group. Therefore for t > 0, macroscopically 
the system is motionless. If we still use ^x^^^(t), u^^^(t)^ for the reverse flow then 
the measure converges weakly to 


(4.1) Mt{d-K,d-v) = { 
with 


^ At (dx) (g) 5(0,1) (dv) + ^ A_t (dx) 0 <5(0-1) (d'v) t <0 


2—V” - / ’ 2 

[ Ao(dx) (g) 5(0,0) (<^v) 


,, . ■ lMdx) + lA_t{dx) t<0 

/Ut(dx) = < 4 2 


(4.2) 


Ao(dx) 


t > 0, 


u(t, x) = 
and decreasing energy; 


XQt(x) • (0,1)+XQ_t(x) • (0,-1) t<0 

0 t > 0, 


(4.3) 


|up/it(dx) = / |vpMt(dx, dv) = 


1 t < 0 
0 t > 0, 


t > 0, 


cf. [BN], Defintion 2.1. 


Remark 4.1. This describes two fronts approaehing eaeh other up until t = 0, 
when they merge and stay at rest, see Figure 7. In the context of the pressureless 
Euler system this is a “sticky” macroscopic solution, cf [BN]. Rather than using 
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particle systems with adhesion dynamics, here we obtain the solution as the limit 
of Hamiltonian dynamics with repulsive foree. We also provide an explanation for 
the loss of energy: all the energy is transferred to a macroscopieally invisible part 
of the system. 


4.2. Transverse flow. It is known that merely requiring decreasing energy does 
not guarantee uniqueness of measure solutions to the system (3.54), see [BN]. This 
persists when comparing the flow of the previous section with the limit of a trivial 
Hamiltonian flow: for this we take the iV-system to consist of molecules that stay 
far enough from each other so that they never interact. We obtain a solution to 
the system (3.53) that coincides with (4.2) for all t < 0. But at t = 0, the moment 
the two fronts meet, instead of merging and staying at rest, they go through each 
other. 

More precisely, for each N = 2n £ N, j = 1,2,N, let 


(4.4) 




/ 

(0,1) if j odd 

< 

(0, —1) if i even. 


For t G M the orbits 


(4.5) 




satisfy the Hamiltonian system (2.1) provided that the interaction range is suf¬ 
ficiently short, for example, a < IjN. (Notice that cjat in Theorem 3.1, and 
therefore in Section 4 satisfies < l/H.) Recalling definition (2.3), set 


(dx,dv) = 




(dx, dv). 


By the definition of Riemann integral, for any continuous bounded /(x, v) we have 


(4.7) 

lim [ /(x, v)M/^^(dx,dv) = i / f {x,t,0,l) dx + ^ [ f {x,-t,0,-1) dx. 

2 Jq 2 Jq 
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1 


t <0 


t = 0 


~r 

t > 0 


Figure 8. Macroscopic flow of (4.9). 


Therefore 


Mt{dx,d'v) 

(4.8) ^ 


: = -At{dx) 0 5(0,1) (dv) + -A_tidx) 0 5(o _i) (d 


V 


The macroscopic density and velocity are 
(4.9) 


= ^^t{dx) + ^A_t{dx), 


u(t,x) = XQti^) ■ (0,1) +Xq_,(x) • (0,-1), t G M, 
see Figure 8. It is easily checked that (3.50), (3.51) hold, and that for all t / 0 


(4.10) 


V (g) vMt,x(dv) = u (g) u. 


Therefore {nt{dx),u(t,x)) also solves weakly the pressureless Euler system for 
t G M. Since / m%idx) = 1 except for t = 0, we can alter 5 at time t = 0 so 
that 


(4.11) 


|u|Vo(dx) = 1, 


still solving equation (3.53). If we still use fit{dx),u{t, x) for the modified solution, 
we then have constant macroscopic kinetic energy in time: 


(4.12) 


|up^i((ix) = 1, t G 


Clearly for all t < 0, (Jlt{dx), u{t, x)), modified or not, coincides with (^((dx), u(t, x)). 
Macroscopically, the same two fronts are approaching each other and, unless we 
know their microscopic origin, we are not be able to tell what will happen for t > 0. 
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Remark 4.2. Notice here the total macroscopie energy of the limit system is 
eonserved in time: 


(4.13) 


|vpMj(dx, dv) = 1, t G M, 


and the macroscopie kinetie energy / \u\^fit{dx) is only part of the total energy 

Jb? 

in general: 

( 4 . 14 ) [ \v\‘^Mtidx,dv) = [ | 5 |Vt(rfx)+ / \v - u\‘^Mtidx,dv). 

JR4 JlR2 J]R4 

Let h{t) = [ |v — u|^ Mt((ix, dv). Then 


(4.15) / |up^t((ix) +/i(t) = 1, t G R. 

4k2 

Notice that h{t) = 0 when t / 0 and /i(0) = 1. Therefore fort < 0, all the energy of 
the system (4.8) is macroscopic kinetic energy which becomes h{0), the fluctuation 
energy, at t = 0. For t > 0 all the energy is again macroscopic kinetic energy. 

S, (4.3). ,4. 

is deereasing in time. Trivially, the eorresponding fluctuation energy h{t) = 0 for 
all t G R. 


Remark 4.3. It is possible that from a Statistical Mechanics point of view the non¬ 
uniqueness described here can he avoided by excluding a set of flows Mt negligible 
with respect to some probability measure. Notwithstanding this, our aim here is to 
understand specific non-uniqueness examples. 

5. Non-Uniqueness from Moments of Measures Satisfying Identical 

Transport Equations 

Section 4 has shown non-uniqueness by comparing moments of the two limit 
flows Mt{dx,d'v) of (4.1) and Mt{dx,d'v) of (4.8). Note that Mt satisfies weakly 
the transport equation 


(5.1) 


dtMt V • VxMt = 0, 
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while Mt satisfies the same with a nonzero kick at t = 0: 

(5.2) dtMt+ = {Mq+- ® 5o{dt), t G M, 

for Mg± = lim Mt. In this section we present two examples where two differ- 
t^>o± 

ent measures solve the same transport equation (5.1), give identical macroscopic 
density and velocity at t = 0, but the macroscopic density and velocity evolve 
differently to provide a non-uniqueness result for the Cauchy problem of the com¬ 
pressible Euler system in space dimension one. 


5.1. Finite systems with velocity exchange. For systems in space dimension 
1, we use identical molecules that move freely until they collide. The arguments 
in this section also hold for systems (2.1) of (finite range, at least) interactions, 
rescaled as in (2.2). In fact, there exist ctat’s such that, for space dimension 1, the 
limit of elastic collisions coincides with the limit of rescaled interactions, see [X]. 
However, such (Tv’s might be too small for the rescaled interaction model to be 
physically better than elastic collisions. For simplicity then, we shall use elastic 
collisions. The complications of finite range interactions were evident in Section 3. 

In the elastic collision model collisions are instantaneous. Momentum and en¬ 
ergy are conserved. Here it will be enough to consider only two kinds of collisions, 
both compatible with finite range interaction dynamics: 


(1) Binary collisions with incoming velocities ui, V 2 and outgoing velocities 
v '2 satisfying 


(5.3) 


Vl+V2 = u'l + V 2 

Vl+Vl = {v[f + {V'2? 


^Vl= v' 2 , V2 = u'l, 


i.e. the molecules exchange velocities (as they are not allowed to go through 
each other). 

(2) Triple collisions, consisting of two molecules exactly as in item (1) and a 
third molecule in between that stays motionless. 




Figure 9. The collisions of subsection 5.1. 


As Zemlyakov shows in his delightful article [Z], several important questions for 
such systems can be answered using the graphs of the molecule positions as func¬ 
tions of time. Following this, the two types of collision we consider are shown in 
Figure 9. 

Consider a 1-dimensional point system , k = 1,..., N obeying 

elastic collision dynamics. Fix any T E (0,oo). For all t G [0,T], assume that all 
collisions are binary or triple as above. 

Proposition 5.1. Let St{x,v) = {x + vt,v). For all t G [0,r] the empirieal 
measures 

1 ^ 

(5.4) Mt{dx, dv) = -Yl 

k=l ' ' 

satisfy 

(5.5) Mj:^\dx,dv) = StMQ^\dx, dv). 

Proof. Merely notice that for each t 

1 iV ^ N 

k=l ' ' k=l ' ' 

since there is a bijection, if multiplicities are taken into account; 


(5,7) {(xr«).T’(*))}«{(T’(o)+«‘r(o).4''> 
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Indeed, the exchange of velocities between the moving molecules of a collision 
establishes a bijection between the orbits before and after that collision. It¬ 
erating this finitely many times brings us back to the initial orbits given by 

The following Lemma will be used repeatedly. 

Lemma 5.2. Suppose that 

(5.8) M^^\dx, dv) = SiMq {dx, dv), MQ^\dx, dv) ^ MQ{dx, dv). 

Then M^^\dx,dv) ^ StMQ{dx,dv). 

Proof. Use the definitions of weak convergence and push forward under St- □ 

As it is standard that Mt{dx,dv) = StMQ{dx, dv) solves weakly the free trans¬ 
port equation 

(5.9) dtMt + vdxMt = 0 
we shell refer to it as the a free transport flow. 

5.2. Euler system from free transport flow. We find here conditions that 
imply that averages with respect to free transport flow satisfy the compressible 
Euler system in dimension 1. The next two subsections provide examples satisfying 
such conditions. 

Lemma 5.3. Suppose that Mt{dx,dv) = StMQ{dx,dv). Then for all (p{t,x) G 
([0,T) X R) and g{v) such that vg{v) G (Mq), we have 

j [ [dtf){t,x)g{v)+dx(j)it,x)vg{v)]Mt{dx,dv)dt 

(5.10) 

-|- / 4>{0,x)g{v)Mo{dx,dv) = 0. 

./r2 


Proof. Straight forward calculation using the definition of the push forward under 
St and the assumption that (f is compactly supported. □ 
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Disintegrating Mt{dx, dv) of Lemma 5.3 as 


( 5 . 11 ) 


Mt{dx,dv) = / Mt^x{dv)iit{dx), 


and for 


( 5 . 12 ) 


g{v){t,x)= / g{v)Mt,x{dv), 


( 5 . 10 ) becomes 


( 5 . 13 ) 


dt(j){t,x) g{v){t,x) + dx<p{t,x)vg{v){t,x) fit{dx)dt 


+ / 4 '{ 0 ,x)g{v){ 0 ,x)noidx) = 0 . 


To apply Lemma 5.3 for g{v) = l^v, and assume G {Mq). Noting that 


( 5 . 14 ) 


u{t,x) = v{t,x) = / vMt^xidv), 


and using the notation 

( 5 . 15 ) ^‘^{t,x)= f {v - u{t,x)fMt^xidv), ^ 3 (t,x)= [ {v - u{t,x)f Mt^x{dv), 

M </ M 

it follows that 


( 5 . 16 ) 


u 2 (t, x) = u^(L x) + ^ 2 (t^ x), 

v^{t, x) = v?{t, x) + 3 u(t, x)^ 2 (L x) + x). 


Then ( 5 . 13 ) for g{v) = l,u, and gives 


( 5 . 17 ) 


{dtcj) + dx 4 >u)iJ,t{dx)dt + / (p{ 0 , x)no{dx) = 0 , 


dt(l)u + dx ( (jiv? + ) ) gt{dx)dt + I (j){ 0 , x)ugo{dx) = 0 , 


1 2 I79 

J_F2 


2 2 


1 2 3^ 

J_F2 


dt 4 > ' 7 ;'^ + 7 ;^^ ]+ dx(l) ' + 17 f d‘t{dx)dt 


2 2 


+ f '('(O,®) Qu^( 0 ,x) + ^0(^3;) = 0 . 
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Moreover, if fit{dx) = p{t,x)dx, ^^{t,x) = 0 and for e{t,x) = 
(5.17) shows that p,u,e solve weakly the Cauchy problem 


C‘^{t,x) 


P = 2/36, 


(5.18) 


dtp + dx{pu) = 0 

dt{pu) + dx {pv?) + dxP = 0 

dt (py + (p^ (y + +7'^^^ = 0, 


P = 2/36, 

p\t=o = p{0, x), u\t=o = u{0, x), e\t=o = e(0, x), 

\ 

the one dimensional Euler system, cf. [CF], p.7. In summary, we have shown: 

Proposition 5.4. For Mt{dx,dv) = StMQ{dx,dv), suppose that G L^{Mq), 
Pt{dx) = p{t,x)dx, and ^^{t,x) = 0. Then p{t,x),u{t,x),e{t,x) as defined above 
is a weak solution to the one dimensional Euler system (5.18). 

The definition of initial conditions for weak solutions here is compatible with 
the one in [dP], p. 2 and [VF], §VII.10. Two examples satisfying the conditions of 
this proposition now follow. 

5.3. Two-layer system. For N fixed, consider N = 2n point molecules xi, 
X 2 ,..., xn on the real line, with 


(5.19) 


, , k , , 1 for A: odd 

a^fc(O) = Y’ = { 

—1 for k even. 


Let the system evolve as in subsection 5.1. After the first n simultaneous collisions 
take place the molecules with labels 1 and N move with velocities 1 and —1, re¬ 
spectively, without ever interacting with any other molecule again. The remaining 
molecules now form a replica of the initial system, reduced by two molecules. 

As in [Z], the graphs of the positions as functions of time show the evolution of 
the system, Figure 10. For 


k=l 


(5.20) 
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X 



Figure 10. Microscopic evolution of subsection 5.3. 
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Figure 11. Macroscopic evolution of subsection 5.3. 
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t 

Figure 12. Microscopic evolution of subsection 5.4. 

5.4. Three-layer system. Consider now for each N = 3n a second system, con¬ 
sisting of N molecules xi,X 2 , ■ ■ ■ ,X]\f on the real line with 

Xk{0) = -^, k = l,...,N, 

for k = 3m — 2 
0 for A: = 3m — 1 

—\/6/2 for A = 3m, m = l,...,n, 

also evolving under elastic collisions as in section 5.1. 

The evolution of the system initialized by (5.28) is shown in Figure 12. Again, 
if for the current system 

~ 1 ^ 

( 5 . 29 ) m'’"’ 

k=l 

by Proposition 5.1, 




(5.30) 
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1 


O' ■ 

t = 0 I 

t = A 

t = .6 

Figure 13. Macroscopic evolution of subsection 5.4. 


On the other hand, as iV ^ oo, 


MQ^\dx,dv) ^ MQ{dx,dv) 

(5.31) 1 1 

= X[o,i]{x)dx (^^-VE/ 2 idv) + -5o{dv) + -d^/^idv) ] , 


therefore 


(5.32) 


^Mt = StMo, N^oo. 


It is again a straightforward calculation that 


(5.33) 


Mt{dx,dv) ^_^^^{x)dx 0 6_^{dv) 

+ ^X[o,i]{x)dx ^ Soidv) + ^^:^^-^{x)dx 0 6 ^{dv) 


Mf describes three layers, each of total mass 1/3, initially overlapping on the 
interval [0,1]. Two of them move with velocities ±\/6/2 for t > 0, while the third 
stays at rest, see Figure 13. The macroscopic density, velocity and energy density 
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given by Mt are 


(5.34) 


p{t,x) = 5^[0,1](*) + 

u{t,x) = -^ 


e{t, x) = 


p{t, x) 


When p{t,x) = 0, take u{t,x),e(t,x) = 0. Notice that 

(5.35) ^3 (0,x)= [ {v-u{t,x)fMt^^{dv) = 0. 

Jr 


By Proposition 5.4, {p, u, e) is also a solution to the Cauchy problem (5.27), clearly 
distinct from the solution {p,u,e). 


Remark 5.5. It is well known that weak solutions to systems like (5.27) are not 
unique, see [D]. This section provides a microscopic interpretation of such macro¬ 
scopic non-uniqueness, showing that such phenomena are quite natural from a 
Hamiltonian point of view. 


Appendix : Motion in a Central Field 

We establish some facts for the motion in dimension 2 of a single particle in an 
external field of potential energy <f> of finite range a: 

(A.l) x"(t) = -ci>'(|x|)^. 

To accommodate (2.2), assume that <1> : (0,oo) —)• [0, oo) satisfies 

(A.2) lim <l>(r) = +oo, <!>' < 0, > 0, 4>(r) ^04^0<r<a. 

1—>0 

Consulting Figure 14, let O be the center of the potential <1>. A molecule m 
enters the range of <k at A with velocity v and leaves at B. For D the middle of 
AB, the path of m in the range of <k is symmetric about OD, by the reversibility of 
the equations of motion. Decompose v{t) into vi(t) and V 2 (t) along AB and OD, 
respectively, and let E be the intersection of OD and the trajectory of m. When 
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a 



Figure 14. Motion in a central field of finite range a. 


m crosses OD it has moved d on the direction of OD. If 9 is the angle between v 
and AB and C is the point on OD with AC of direction v, then 

(A.3) d = DE < CD = AC ■ sin0 < AO ■ sin0 = iTsin0. 


Let T be the time it takes m to travel from A to B. 

4(7 

Lemma A.l. For a the range o/<h, T and v as above satisfy T < —. 

V 

Proof. From (A.l), 


(A.4) 




X • X 


X • 


Xn 


X • X 


X2 - $'(|x|)^ + ^'(|X|j^-X2. 

X xL 


For X 2 < 0 and as —;— < 0 for f G (0,T/2), and as is convex, the first term of 
dt 

this is negative and, if xi is also negative, the sum of the remaining two terms is 
also negative provided that 


(A. 5 ) 


—x'2\yd\^ + X2(x • x') > 0 —x'2Xl + X2x'i < 0, 
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since is decreasing. Now note that —X 2 X 1 + X 2 X 1 stays constant in time and the 
inequality is satified at t = 0. Therefore V 2 is concave and by (A.3) 

(A.6) ■ — < d < asinO, 

2 2 

which, along with U 2 ( 0 ) = usin0, concludes the proof. □ 

Still in Figure 14, let ZACD = (j). Denoting the distance of O from AC (the 
impact parameter) by a, by [LL], p. 49^ 


(A.7) 


</>(«) = [ 
Jr, 


a 


rmn 2 . I 1 ^ 


dr, 


Ht/1- —- „ 


where E = -mv"^ and rmin is a zero of the radicand: 


(A.8) 


1 - 




E 


= 0 . 


Lemma A.2. For interaction potential as in (2.2), rmin = rminict) is increasing 
and 4'{a) is eontinuous on [0,oo). 


Proof. For fixed a and E, the function 
(A.9) 


o? <h(r) 
E 


is strictly decreasing from +oo to 0 for r > 0 and the pre-image rmin of 1 satisfies 
(A.8), or 


(A.IO) 


,, ^{Tmvn)\^'’^ 

a = 1 1--— 1 

tj 


showing that a = airmin)-, and therefore rmin = rmin{oi), is increasing. 


^Note here that [LL]’s analysis of motion in a central field in their §14 is valid for any central 
field, including the ones with finite range. 
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To show that (j) is continuous, change the variable in (A.7) via r = rminU' 

1 


(A.ll) 


(j){a) = 


by (A. 10) 


y2 ^ l_ 




^{rminy)\ 1 


zdy 


1 


'1 E-^{rminy) 1 


E ) y- 
--dy. 


E - <l)(r„ 


y^ 


From (A.8) we have 


(A. 12) E > 

and since <h(r) is decreasing, 


(A.13) 




therefore 

(A.14) 

with 

(A.15) 


1 


IE ^(j'miny') 


< 


1 


E - ^{rmin) y"^ 


1-- 

y 



TT 

2‘ 


In other words, the integrand of cj) is dominated by an integrable function. This, 
and the continuity of r^in in ci; show that (j) is continuous in a. □ 


Corollary A. 3. For any 0 < 4>o < 7r/2, there exists 0 < ao < a such that 
0(«o) = 00- 

TT 

Proof. Just use continuity and that i;A(0) = 0 (“head-on collision”), (j){cr) = — (no 
interaction). □ 

As is well known, motion in a central field also describes a system of two bodies 
interacting with each other via <I>, a function of their distance, in a coordinate 
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system with its origin at the center of mass of the system. The formulas for this 
transformation are in [LL], §13. 
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